Simulation of relativistically colliding laser-generated electron flows 
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The plasma dynamics resulting from the simultaneous impact, of two equal, ultra-intense laser 
pulses, in two spatially separated spots, onto a dense target is studied via particle-in-cell (PIC) 
simulations. The simulations show that electrons accelerated to relativistic speeds, cross the tar- 
get and exit at its rear surface. Most energetic electrons are bound to the rear surface by the 
ambipolar electric field and expand along it. Their current is closed by a return current in the 
target, and this current configuration generates strong surface magnetic fields. The two electron 
sheaths collide at the midplane between the laser impact points. The magnetic repulsion between 
the counter-streaming electron beams separates them along the surface normal direction, before 
they can thermalize through other beam instabilities. This magnetic repulsion is also the driving 
mechanism for the beam-Weibel (filamentation) instability, which is thought to be responsible for 
magnetic field growth close to the internal shocks of gamma-ray burst (GRB) jets. The relative 
strength of this repulsion compared to the competing electrostatic interactions, which is evidenced 
by the simulations, suggests that the filamentation instability can be examined in an experimental 
setting. 

PACS numbers: 52.38Kd, 52.38.Fz, 52.65.Rr 



I. INTRODUCTION 

The dynamics of the transient electromagnetic field 
driven by relativistic electronic surface currents, which 
are generated by the interaction of ultra-intense laser 
beams with solids [J_H3] , are of relevance to important ap- 
plications such as laser-driven ion acceleration [U [5] and 
fast ignition [6]. The surface magnetic and electrostatic 
fields can be of the order of - 10 3 T and - 10 12 V/m [7HS], 
respectively, during the ultra-intense (/ > 10 19 W/cm 2 ) 
laser-solid interaction, and such fields can last several 
tens of picoseconds after the laser pulse ends. Charging of 
a laser-irradiated target, which is attributed to the escape 
of hot electrons generated during the laser-plasma inter- 
action, has been detected in several experiments [TO] [11] 
using proton probing techniques. Because of their con- 
finement by the sheath electric field and surface magnetic 
field, most of the hot electrons are bound to the target 
surface (3] IT21 - H4] as they move away from the interac- 
tion region. The surface current is closed by the return 
current within the target and this current loop encloses 
a strong magnetic field. This magnetic field generation 
mechanism is known as the fountain effect. Recent ex- 
periments [151 116) have demonstrated that, during the 
ultra- intense laser irradiating the target, the magnetic 
field tied to the electron current expands along the solid's 
surface at a speed close to that of light. 

If two ultra-intense laser pulses of comparable intensity 
irradiate the target's surface at two separate points si- 
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multaneously, they produce two relativistic electron flows 
at the rear surface. Both flows collide head-on half-way 
between the two interaction points at the planar rear 
surface. The superposition of their respective currents, 
which are oppositely directed and comparable in mag- 
nitude, initially results in a net current density, which 
is low compared to the current density of each electron 
beam. This reduction of the current density results in 
a weakened magnetic field in the beam overlap region. 
The time-evolution of the counter-streaming beams and 
of the electric and magnetic fields, which result in their 
thermalization, depends on the nature of the instabilities 
and processes that are involved. 

Previous related experiments have examined the simul- 
taneous interaction of two or more laser pulses with a 
solid target. The laser pulses had intensities of 10 14 — 
10 15 W/cm 2 [T7J [T5J. In these experiments, the expan- 
sion of plasma bubbles at the front surface of the target 
was examined with the aim to study magnetic reconnec- 
tion of the megagauss-fields, which encircle the individual 
plasmas and confine them. The moderate laser intensity 
has limited the particle flow speeds to the nonrelativistic 
regime. Such flows tend to thermalize electrostatically. 

Here the speeds of the electron sheaths, which form 
before ion jets are launched, are relativistic and the mag- 
netic fields sustained by their currents are confined to a 
narrow layer close to the surface. The instabilities that 
result in the thermalisation of the relativistic fast electron 
flows differ from the nonrelativistic ones in Ref . [T7] [TS] 
that involve electrons and ions. 

The counter-streaming relativistic electrons give rise to 
a wide range of beam instabilities in this overlap layer. 
The potentially large number of instabilities can be re- 
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duced to three, if we assume that only the electrons in- 
teract, that the thermal spread of the electron speeds 
is smaller than the difference of the mean speeds of both 
beams and that magnetic field effects in the beam overlap 
region are negligible. Negligible here means that the elec- 
tron gyrofrequency is small compared to the growth rates 
of the instabilities. We obtain under these conditions the 
two-stream instability, the oblique mode instability and 
the filamentation or beam-Weibel instability. These in- 
stabilities are discussed in detail in Ref. 19 . Here we 
summarize their properties. 

Let k and be the wave vector of the growing waves 
and the difference vector between the mean velocities of 
the counter-streaming electron beams, respectively. Let 
a be the angle between k and v^. The two-stream insta- 
bility [20] is purely electrostatic for k || v;, (a = 0). ft 
saturates by the formation of electron phase space holes 
in the nonrelativistic [21] and relativistic regimes [19]. 
The coalescence instability [2T] triggers the collapse of 
these nonlinear structures, which thermalizes the elec- 
tron flow. What remains from the two counter-streaming 
electron beams is a single hot electron distribution with 
a low mean speed modulus. The net current is small and 
the magnetic field it drives is weak. The two-stream in- 
stability grows slower than the oblique mode instability 
if |v;,| ~ c. The oblique modes have a mixed polarity and 
they are almost electrostatic. Their wave vectors k form 
angles < a < 90° with respect to v?, and they become 
two-stream unstable modes if a = 0. They also saturate 
by electron trapping [22]. Both instabilities thermalizc 
the electrons electrostatically. 

If the counter-streaming electron beams have a similar 
density and temperature and if they collide at a moder- 
ately relativistic speed, then both electrostatic instabil- 
ities are outgrown by the filamentation instability |23j . 
The latter is driven by the magnetic repulsion of elec- 
trons, which move into opposite directions. It saturates 
by magnetic trapping [24) and its final state are current 
channels, which contain charged particles that have the 
same current direction. These channels have a diameter 
of the order of the electron skin depth and are surrounded 
by strong magnetic fields. The plasma flow is eventually 
thermalized through the repeated mergers of the current 
channels. However, the magnetic fields remain strong in 
a broad spatial interval [25H25] . 

Our 2D PIC simulation, which uses the 2D3V code 
LAPINE [2J5], resolves a cross-section of the target that 
contains the surface normal. The electron flow at the 
rear surface is thus confined to a layer close to the one- 
dimensional surface with a width that is comparable to 
an electron skin depth. The colliding electron sheaths 
can drive two instabilities in this geometry. The two- 
stream instability, which results in electrostatic oscilla- 
tions along the surface direction, and the filamentation 
instability, which separates both electron beams along 
the surface normal. Our simulation demonstrates that 
for the selected initial conditions the magnetic interac- 
tion is more important than the electrostatic one. A 



laser-plasma experiment with similar initial conditions 
would thus allow us to drive the filamentation instability 
on the two-dimensional surface of the target. The experi- 
mentally measured data can be compared to the existing 
2D and 3D simulation studies (See for example [251 - 125] ') of 
colliding leptonic flows and test their validity. Studying 
this instability with a controlled laboratory experiment 
would provide us with a much needed better understand- 
ing of the magnetic field generation within gamma-ray 
burst (GRB) jets, which are spectacular releases of ener- 
getic electromagnetic radiation at cosmological distances 
[50] . A direct comparison of the filaments observed in a 
controlled laboratory experiment and those believed to 
exist close to the internal shocks of GRBs is possible if 
the plasma is collisionless, because then all observables 
can be scaled between both systems with the help of the 
electron plasma frequency |31) . 

This paper is structured as follows. Section 2 discusses 
the numerical scheme and the initial conditions. The 
results are presented in Section 3 and they are discussed 
in Section 4. 



II. SIMULATION MODEL 

Particle-in-cell simulation codes solve the numerical 
approximation of Maxwell's equations on a grid and they 
approximate the plasma by an ensemble of computational 
particles. The particle momenta are updated with the 
help of the relativistic Lorentz force equation. This sim- 
ulation method is discussed in detail in Ref. [32]. We 
restrict our simulation to a 2D geometry. 

Our initial and boundary conditions are as follows: 
Two equal laser pulses hit a thick plasma slab simul- 
taneously and at normal incidence as depicted in Fig. [I] 
The target with a width of HOAo and thickness of 5Ao 
is initially located between z = 10Ao and 15Ao, where 
Ao = ljUm is the laser wavelength. The target consists of 
an initially charge-neutral mixture of electrons and Al 10+ 
with mass — 27m p , where m p = 1836m e is the pro- 
ton mass. The initial density of the target is set to 50n c , 
where n c = 1.12 x 10 21 cm~ 3 is the critical density. The 
initial temperatures of the electrons and ions are both set 
to 1 keV for computational reasons. The simulation box 
size L z x L y — 30A x 110A is resolved by 1500 x 5500 
cells, and the electrons and Aluminium ions are repre- 
sented by 64 and 16 computational particles per cell, re- 
spectively. Two p-polarized laser pulses are incident nor- 
mally from the left boundary and focused on the target 
at y = ±30 Ao, respectively. The midpoint between the 
laser pulses defines y = 0. The pulses rise up in the first 
4To with a Gaussian profile and then maintain the peak 
intensity for 40T , where T ~ 3.3 fs is the laser cycle. 
The maximum intensity of the pulses is 5 x 10 19 W/cm 2 , 
corresponding to ao = 6. The laser pulses have a Gaus- 
sian spatial profile with a spot radius of 5Ao- The time 
step is 0.007T . For both the transverse and longitudinal 
boundaries, absorbing boundary conditions are used for 



3 



Laser 1 




11M 



Laser 2 



y a 



5 ^0 



FIG. 1: (Color online) The simulation setup: The target 
(blue) has a thickness of 5Ao along the z-direction and is lo- 
cated between z — lOAo and 15Ao in the simulation box. The 
target is hit by two equal laser pulses (red). 




14 16 18 20 22 24 



14 16 18 20 22 24 



the fields and particles. 



III. RESULTS 

The impact of the double pulse generates two rarefac- 
tion waves, as can be seen in Fig. [2j The electrons are 
accelerated mainly by J x B heating (i.e., by the os- 
cillating component of the laser's ponderomotive force) 
[33] to a temperature of ~1.0 MeV as the laser pulses 
irradiate the target, which is close to Haines' scaling 

T h = m e c 2 (\/l + V2a - 1) = l.OlMeV [34j with a ve- 
locity of 0.95c. The electron temperature obtained here 
is much lower than that given by ponderomotive scaling 
(~2.6 MeV) [35]. This discrepancy can be attributed to 
the fact that the electrons only interact with the laser 
pulse during a fraction of the laser period before being 
accelerated forward beyond the laser penetration region, 
due to that in our simulation that ultra-intense laser 
pulses irradiate an overdense target with a steep density 
interface. This situation is different from the mechanisms 
proposed by Kemp et al. [3S] and May et al. [37]; a l° w 
density shelf is required for the former and the electron 
energy is gained in the vacuum from the transverse laser 
field in the latter case. 

The hot electrons propagate through the target and 
form the two high energy density channels visible in Fig. 
[2jb) at y = ±30 Ao A small part of the energetic elec- 
trons escape into the vacuum at the rear surface of the 
target. Their current, visualized in Fig. [2ja), is not bal- 
anced by an ion current and generates a strong sheath 
electromagnetic field, which reflects most of the hot elec- 
trons and confines them. The sheath field accelerates ions 
on the target's surface (See Fig. [2^c)), whose larger iner- 
tia implies that they trail the electrons in Fig. J2jd) . The 
plasma density of both parabolic expanding rarefaction 
waves decreases with increasing distance from the target 



FIG. 2: (Color online) The amplitude of the electron current 
density (j J y +i* J z \) in the simulation plane (a), the electron 
kinetic energy density in units of m e c 2 n c (b), the ion charge 
density in units of en c (c) and the electron charge density in 
units of — en c (d). The amplitude of the current density is 
in units of en c c. All color scales are 10-logarithmic and the 
simulation time is t=120To. 



and the charge separation drives an ambipolar electric 
field. 

The flow of hot electrons at the surface of the target 
and the return current, which flows within the target to 
provide current closure, result in the growth of magnetic 
fields orthogonal to the simulation plane, which is here 
the x-direction. In the case of a single laser pulse, strong 
surface magnetic fields will expand from the exit point of 
the electrons at the rear end of the target until they cover 
the surface uniformly. Here the currents driven by both 
laser pulses will eventually collide and a more complex 
magnetic topology is revealed by Fig. [3] The distribution 
of B x at the early time 60T in Fig. Ma), which is trig- 
gered by each laser pulse, equals that observed for single 
laser pulses. The hot electron's surface current sheaths 
expand at the speed 0.7c, but they have not yet collided 
at y = 0. The polarity of B x at y — ±30A switches, 
because the electrons above the interaction point flow in 
the opposite direction compared with those below that 
point. 

The hot electron sheaths have collided at the time 
102To in Fig. [3jb) . The surface magnetic field has de- 
creased in amplitude compared to that at 60T and the 
field is more confined along z. A sharp transition between 
the magnetic fields occurs at y = and a localized mag- 
netic structure is observed at z rj 15.5Ao and y « — 4Ao. 
The polarity (negative amplitude) suggests that it is tied 
to the electron sheath generated by the upper laser pulse. 
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FIG. 3: (Color online.) Transverse magnetic field (B x ) dis- 
tribution at t=60T o (a), 102T (b), 120T (c), and 140T (d), 
respectively. The amplitude is expressed in units of m e coJo/e. 



The current distribution is not perfectly symmetric at 
this time as we may expect from our symmetric simula- 
tion setup. However, the approximation of the plasma 
by computational particles introduces some randomness. 
Even a small difference of the current or charge density of 
the surface electrons can result in significant differences 
in the beam evolution, because their separation by the 
repulsive magnetic interaction corresponds to a plasma 
instability. Statistical variations could explain the ob- 
served asymmetric current distribution. Such differences 
do not limit the comparison of our PIC simulation results 
with that in experiments, since two identical laser pulses 
would never be obtained in experiments. 

This current structure and the magnetic field it drives 
have expanded further in Fig. [3^c) (t — I20T ). A second 
similar one has developed at y > and z > 15Ao- We 
call these structures magnetic wings. The distribution 
of B x at the time 140To in Fig. [3jd) suggests that the 
currents responsible for the wings are about to reconnect 
to the surface currents at \y\ < 30Ao and z w 15Aq. We 
infer this from the magnetic field polarity. We find in 
Fig. [3][d) a magnetic field patch with a positive ampli- 
tude in the interval close to y ~ 40Ao and z w 18Ao and 
in an interval close to y ss 30Ao and z 20 Ao- Both are 
about to merge. Such a merger implies that the currents 
encircling both structures reconnect. An unambiguous 
demonstration of such a reconnection can only be pro- 
vided by the electronic currents. They are, however, too 
noisy for this purpose. The magnetic field distribution 
is smoother, because it is connected to the well-resolved 
total current and not to the fluctuating current densities 
within individual simulation grid cells. 



The magnetic fields observed in Fig. [3] are transient. 
The current, which has been driven directly by the laser 
pulses, will eventually be dissipated. The thermoelec- 
tric instability, an instability driven by a plasma density 
gradient that is not parallel to the temperature gradient 
|38j , will eventually develop and remagnetize the plasma. 
The dynamics of the magnetic fields it drives and their 
reconnection have been examined in the experiments per- 
formed in Refs. [T71 |T5] . We do not consider them here. 

Figure [4] displays J y , which corresponds to the surface 
current along the initial target boundary, and the electric 
field's Eg component at two times (only the target's rear 
side is shown). The ambipolar E z is tied to the density 
gradient of the rarefaction wave and it is here practi- 
cally electrostatic. Strong surface currents are present 
in Fig. Qa) at z ~ 15Ao- The fluctuations observed 
for z < 15Arj are thermal noise. No such fluctuations 
are seen for z > 15Ao, since the lower plasma density im- 
plies lower fluctuation amplitudes. The current filaments 
at z ss 15Ao have a thickness, which is a few times the 
electron skin depth c/uj p w 0.023Ao- Consider the point 
y = 30 A in Fig. |4]Ja). The surface current above this 
point is negative, which implies that electrons move to 
increasing values of y. It has the opposite sign below. 
The hot surface electrons thus move away as expected 
from the point y — 30Ao . In addition, a lower return cur- 
rent density at z < 15Ao is visible. Both currents encircle 
the surface magnetic field visible in Fig. [3jb) . 

The surface current in Fig. [4[a) maintains a practi- 
cally constant strength from \y\ = 30 A up to y = 0. 
The electron flow has not been thermalized by an elec- 
trostatic instability, as this would demagnetize the flow 
around y = and increase the electron temperature. The 
currents have been separated instead along the normal 
direction of the target's rear surface in a small y-interval 
close to y « 0. The surface current driven by the up- 
per laser pulse is magnetically expelled at y m —4\ Q and 
z « 15.5Ao and J y reaches here a value 5.4 x 10 1:L A/cm 2 . 
The total current density is actually stronger, because 
the expelled current has J y and J z components due to 
its deflection along z. Its spatial correlation with the 
magnetic structure in Fig. |3jb) indicates that the ex- 
pelled current drives it. The surface current driven by 
the lower laser pulse initially reconnects with the return 
current driven by the upper laser pulse. We attribute 
again the different behaviour of the electron flows driven 
by the upper and lower laser pulses to the break of sym- 
metry by the plasma approximation by a finite number of 
computational particles and the resulting statistical vari- 
ations. The symmetric wings in Fig. [3^c) evidence that 
both current sheaths get expelled at a later time. Note 
that we can not clearly associate a current with the over- 
plotted contour lines of the wings in Fig. |4]jb), because 
the current is too weak and noisy. The total current is, 
however, preserved since Fig. |3jc) clearly demonstrates 
that the peak amplitudes of the surface current and wing 
magnetic field are practically identical. 

The reason for the widening of the expelled current 
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FIG. 4: (Color online.) The transverse current density (J y ) 
distribution at t=102Tb is shown in (a) and that at 120Tb in 
(b). Several contour lines of B x are overplotted in (b) using 
the same color scale as for the current density. The electric 
field (E z ) distribution at t=90Tb is shown in (c) and that 
at 120Tb in (d). The color scale of the current density is 
restricted within [-0.3 0.3] in order to distinguish the current 
distribution clearly in the vacuum. The current density and 
the electric field amplitude are expressed in units en c c and 
m e aoo/e, respectively. 



sheath, which is responsible for the spreading out of the 
magnetic wing with increasing distance from the target 
surface in Fig. [3^c) , is given by Figs, gc) andgd). The 
parabolic ambipolar electric field visible in both plots 
outlines the front of the two rarefaction waves, which are 
driven by the laser pulses. This interpretation is sup- 
ported by the electron and ion distributions in Figs, gc) 
and gd) . Note that no ambipolar electric field can be 
observed at the original target boundary z = 15A , be- 
cause the strong surface magnetic field suppresses the 
electron mobility along the surface normal. A compari- 
son between Fig. gc) and Fig. gd) demonstrates that 
the magnetic wings are located within the rarefaction 
wave. A current closure is thus achieved by the rarefac- 
tion wave's plasma. The plasma density within the rar- 
efaction wave is below that of the solid target, which 
reduces the supported current density. The currents and 
the associated magnetic fields spread out in space. 

It is interesting to know how the laser power and, thus, 
the characteristic speed of the electron sheath affect the 
interaction of the electrons and the magnetic fields in 
the sheaths' overlap layer. More specifically, we want to 
know if the filamentation instability could also be exam- 
ined with weaker laser pulse intensities, like those used 
in the previous experiments |17l 118] . We have performed 
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FIG. 5: (Color online.) Transverse magnetic field (B x ) distri- 
bution at t=60T (a), 100T (b), 120T (c), and 140T (d), re- 
spectively, for the case with a laser intensity of 5x 10 18 W/cm 2 . 
The other parameters are the same as that in Fig. [3] 



for this purpose a scries of simulations with lower laser 
intensities, i.e., J = 5x 10 16 W/cm 2 , 5 x 10 17 W/cm 2 , 
and 5 x 10 18 W/cm 2 , respectively. The magnetic field 
repulsion in the sheath overlap layer is not observed in 
any of the simulations. Figure [5] shows the evolution 
of the transverse magnetic field for the case study that 
employed a laser intensity 5 x 10 18 W/cm 2 . The mag- 
netic field amplitude driven by the surface current is 
much weaker than that observed for a laser intensity 
5 x 10 19 W/cm 2 (Fig. g. The maximum magnetic field 
here is ~ 2.8 x 10 7 Gauss at t = 100T , only one third of 
the latter (~ 7.4 x 10 7 Gauss). Both the current density 
and the expansion speed of the electron current along 
the surface here (~ 0.4c) are much smaller than that of 
the latter. Most importantly, Fig. [5] reveals a different 
magnetic topology close to y ps compared to Fig. [3] 
The electrons from both sheaths do not interpenetrate 
in the simulation that uses the weaker laser pulse. This 
implies that the driving currents are closed by the re- 
turn current within the target. A magnetic distribution 
develops at the surface, which is stationary during the 
resolved time interval. The absence of counter-streaming 
electron beams implies that no filamentation instabil- 
ity can develop here. If the electron beams would in- 
terpenetrate close to y « 0, their nonrelativistic speeds 
f=a 0.4c observed in the simulation would probably result 
in a thermalization through the quasi-electrostatic two- 
stream and oblique modes |23j . 
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IV. CONCLUSION 

In summary, we have examined the plasma processes 
triggered by the impact of a laser double-pulse on a solid 
target, using PIC simulations. The large 2D simulation 
box has resolved a cross section of the target. Each laser 
pulse has accelerated the electrons to relativistic speeds. 
The electrons have mainly been accelerated along the 
propagation direction of the laser pulse and they have 
crossed the target. The electrons re-emerged on the tar- 
get's rear surface, where they have been deflected by the 
self-generated ambipolar electrostatic field. They have 
expanded along the target's surface and they have col- 
lided at the midpoint between the laser axes. The plasma 
did not thermalize electrostatically when the electron 
sheaths collided. Their relativistic speeds resulted in 
a dominant magnetic rather than electrostatic interac- 
tion. Our 2D geometry implied that the current sheaths 
were expelled from the target by their mutual interaction 
and moved into the rarefaction wave, provided that the 
driving laser pulse is ultra-intense. The resolution of a 
second dimension of the target's surface by a currently 
too expensive 3D PIC simulation would allow for stud- 



ies of the filamentation instability [THl [551 ISO] within the 
overlap layer of the colliding electron sheaths. Our sim- 
ulation results thus suggest that this instability, which is 
thought to magnetize energetic astrophysical flows, can 
be observed in a laboratory experiment provided that the 
laser pulse is ultra-intense and the electron flow speeds 
are relativistic. 
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